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A "drivebelt" stadium billiard with boundary consisting of circular arcs of differing radius con- 
nected by their common tangents shares many properties with the conventional "straight" stadium, 
including hyperbolicity and mixing, as well as intermittency due to marginally unstable periodic 
orbits (MUPOs). Interestingly, the roles of the straight and curved sides are reversed. Here we 
discuss intermittent properties of the chaotic trajectories from the point of view of escape through 
a hole in the billiard, giving the exact leading order coefficient lim^oo tP(t) of the survival prob- 
ability P{t) which is algebraic for fixed hole size. However, in the natural scaling limit of small 
hole size inversely proportional to time, the decay remains exponential. The big distinction between 
the straight and drivebelt stadia is that in the drivebelt case there are multiple families of MUPOs 
leading to qualitatively new effects. A further difference is that most marginal periodic orbits in this 
system are oblique to the boundary, thus permitting applications that utilise total internal reflection 
such as microlasers. 



Intermittency is a general dynamical phe- 
nomenon in which a single chaotic trajectory 
is interrupted at apparently random intervals 
by quasi-regular behavior of varying and un- 
bounded duration. While first described in the 
context of turbulence, it is ubiquitous in nature 
since one of the sources of intermittency, ellip- 
tic islands approached more and more closely 
by the long quasi-regular intervals, are generic 
features of Hamiltonian systems. Another, and 
more mathematically tractable, source of in- 
termittency is that of marginally unstable pe- 
riodic orbits (MUPOs). In contrast to ellip- 
tic islands, MUPOs can exist in systems with 
strong ergodic properties such as the celebrated 
stadium billiard introduced by Bunimovich in 
1974. This billiard consists of a particle con- 
fined to a region bordered by two straight par- 
allel lines and two semicircles tangent to these, 
in which the MUPOs are "bouncing ball" or- 
bits reflecting between the parallel sides indef- 
initely. Orbits with many consecutive bounces 
on the straight segments are quasi-regular, and 
when reaching a curved end of the billiard, they 
return to make a quasi-regular sequence of simi- 
lar length; only after many such sequences does 
the orbit return to the chaotic region of phase 
space. In the case of two circular arcs con- 
nected by non-parallel tangents, the drivebelt 
stadium, orbits between the straight segments 
are of bounded length, but orbits close to pe- 
riodic orbits in the larger arc are now quasi- 
regular. After slowly precessing and reaching 
a straight segment, an analogous return takes 
place in a sequence of similar length in the re- 
verse direction. Now, however, there are po- 
tentially many such families of MUPOs, leading 
to the richer phenomena of multiple intermit- 
tency. Here we quantify this multiple intermit- 
tency by drilling a hole in the billiard and wait- 
ing for the particle to escape. The probability 
of surviving for a specified time, P(t) decays al- 



gebraically, with a coefficient expressed exactly 
as a sum of contributions from the MUPOs; nu- 
merical simulations confirm these calculations. 
Properties of the straight stadium are retained, 
but new possibilities arise from the multiple 
intermittency, which resembles more closely a 
generic mixed system with many elliptic is- 
lands. Applications in atom, electron and op- 
tical billiards are discussed; unlike the straight 
stadium, the oblique reflections of typical MU- 
POs in the drivebelt now allow total internal 
reflection, which allows trapping at the very 
small scales needed for microlasers. These wave 
billiards also raise intriguing questions for the 
future - how is multiple intermittency manifest 
in the quantum regime? 



I. INTRODUCTION 

The study of mathematical billiards [1] has received 
much attention in recent years as it facilitates a bet- 
ter understanding of the underlying dynamics of many 
physical systems while also providing useful mod- 
els with applications in areas ranging from nonlinear 
dynamics and statistical physics to quantum optics. 
Classically, the billiard's phase space is defined by the 
position and momentum variables of a point particle 
moving with constant speed in a straight line inside 
a compact domain Q £ M. d experiencing mirror-like 
reflections at the billiard boundary dQ. As the par- 
ticle trajectories do not depend on the overall energy, 
the billiard dynamics (integrable, mixed or chaotic) 
is completely controlled by the boundary geometry. 
The addition of a "hole" or "leak" in phase space 
through which an initial ensemble of particles may 
escape through was first suggested by Pianigiani and 
Yorke [2 j. This made the survival probability P(t) (a 
monotonically decreasing function of time t) an impor- 
tant statistical observable, able to describe and clas- 
sify the internal dynamics of the billiard or other dy- 
namical systems. Moreover, it motivates the following 
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question: How do the long-time escape properties de- 
pend on the dynamics, the size and position(s) of the 
hole(s)? 

To a first approximation, the exponential rate of es- 
cape in strongly chaotic dynamical systems is propor- 
tional to the size of the hole, inversely proportional 
to the size of the accessible part of the phase space 
and independent of the hole's position. We shall be 
concerned with two dimensional billiards where this 
statement can be simply expressed as 

In P(t) e 
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such that is the mean free path between collisions, 

e is the hole size, and \Q\ and \dQ\ are the area and 
perimeter of the billiard respectively. Equation (1) is 
a good approximation for small holes and times longer 
than the scale set by correlation decay which is pro- 
portional to one over the Lyapunov exponent of the 
chaotic saddle. Exact formulas for the subsequent cor- 
rection terms were derived recently in [3] as sums of 
correlation functions depending on both hole size and 
position, and were numerically tested in the dispers- 
ing diamond billiard with a hole placed in coordinate 
space along the billiard boundary. Holes in momen- 
tum space have also been considered and are of inter- 
est when studying for example emission from dielectric 
micro-cavities 0], The full escape properties of open 
systems are in theory described by the properties of 
the invariant chaotic saddle(s) along which a condi- 
tionally invariant measure concentrates in time [5] but 
are however difficult to compute in practice [6]. Nev- 
ertheless, although placing holes through which parti- 
cles may escape was originally a perturbative approach 
for generating transient chaos [21 E] , it has also been 
successful as a kind of chaotic spectroscopy [5j used 
for retrieving information about the internal dynam- 
ics go El. 

Generic dynamical systems however such as the 
standard map |10j are weakly chaotic and have a mixed 
phase space where Kolmogorov-Arnold-Moser (KAM) 
hierarchical islands coexist with a chaotic sea. Al- 
though the dynamics in the hyperbolic regions of the 
phase space is exponentially mixing, trajectories ap- 
proaching the islands 'stick' there for arbitrarily long 
periods of time before they are re-injected back into 
the 'deep' hyperbolic regions of the chaotic sea. The 
sticky dynamics [TT] described above is closely associ- 
ated with the more general phenomenon of intermit- 
tency and is captured by a cross-over from exponen- 
tial to an asymptotic power-law decay in the survival 
probability P(t) ~ t~ a , with a > 1, and is indepen- 
dent of the hole's position, assuming that it is placed 
well inside the chaotic part of the phase space. In 
recent studies [12], this was interpreted as an effec- 
tive splitting of the chaotic saddle into hyperbolic and 
nonhyperbolic components. The question of a univer- 
sal power-law exponent in mixed phase space systems 
is still an open problem though recent results suggest 
a w 0.57 [T3]. 



Power-law asymptotic decays are not restricted to 
mixed phase space systems. For example, power-law 
tails of order ~ i _1 were observed in the famous Buni- 
movich stadium billiard [14] . a paradigmatic example 
of chaotic billiards. The stadium billiard [IS] is com- 
posed of two semicircles joined smoothly (C 1 ) by two 
parallel straight lines. The defocusing mechanism of 
chaos present guarantees a positive Lyapunov expo- 
nent (exponential separation rate of nearby trajecto- 
ries) almost everywhere, the exception being a zero- 
measure continuous family of period two marginally 
unstable periodic orbits (MUPOs) called "bouncing 
ball" orbits trapped forever between the stadium's par- 
allel walls. Bouncing ball orbits have been shown to 
lead to an intermittent, quasi- regular behavior which 
effectively causes the closed stadium to display some 
weaker chaotic properties such as an algebraic decay of 
correlations |16j and a non-standard central limit the- 
orem (CLT) requiring a tlogf normalization [T71 118j . 
Note that "whispering gallery" orbits which slide al- 
most tangentially along the curved boundaries are only 
intermittent in discrete time but not in continuous 
time as their maximum length is bounded. Quantum 
mechanically bouncing ball orbits are known to cause 
scarring |19) . deviations from Random Matrix Theory 
(RMT) (especially in the A 3 -statistics) if not treated 
appropriately [201 [21] , while the system is not quan- 
tum uniquely ergodic [22) and an h dependent 'island 
of stability' appears to surround them [23]. Therefore 
MUPOs are important in understanding the quantum 
to classical correspondence [23] of the phenomenon of 
stickiness. 

It was recently found by the authors that by po- 
sitioning a hole on the boundary as to intersect the 
continuous family of bouncing ball orbits it is possible 
to explicitly single out long surviving orbits hence cal- 
culating the constant C in closed form as a function of 
hole position [25] hence supporting the interpretation 
of Ref.[T2]. Unlike past investigations of open stadia 
[2"fj] , this result included the exact coefficient of the sur- 
vival probability, allowing for the accurate prediction 
and hence optimization of escape time distributions 
while also providing answers to the aforementioned 
motivational question. In this paper we review and 
generalize the approach of [25] to circular type MU- 
POs in the context of the drivebelt billiard [T51 12TH2T)] . 

The drivebelt billiard is constructed by two circles of 
radii r and R > r with centers displaced by a distance 
d. The circles are then connected by their common 
outer tangent lines of length L (see Figure [l]) and the 
interior arcs are removed such that the boundary is 
C 1 smooth and convex and d = y/L 2 + (R-r) 2 . For 
the purpose of our current investigation, only half the 
arc length of the larger circle <f> £ (vr/2, tt) needs to be 
known as shown in Figure [T] Note that in the limits 
<j) — > 7r/2 and <f> — > n the stadium billiard and circle 
billiards are recovered respectively. 

This construction, also known as a 'tilted' stadium, 
or a 'squash' was originally proposed by Bunimovich 
(unpublished), investigated numerically and experi- 
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mentally using ultra cold atoms confined by lasers 
beams in [27], and later rediscovered in [16] and [28] . 
where it was first called a 'drivebelt' billiard and its 
polynomial mixing rates were studied in detail. Re- 
cently, Balint et al. have studied a two-parameter set 
of two-dimensional billiards with one of the limiting 
cases being the drivebelt and conjecture that they ob- 
tain ergodic non-dispersing billiards which are close to 
their drivebelt limit [25]. The billiard is hyperbolic, er- 
godic, and Bernoulli |15j and remains chaotic no mat- 
ter how short L is, but turns into the integrable cir- 
cle when L — > 0. The mechanism of chaos present is 
the defocusing one (same as for the stadium billiard) 
where a parallel beam of rays becomes convergent and 
then divergent after traveling twice the focal length. 
Therefore expansion in phase space is guaranteed al- 
most everywhere |30j . the exception here being a set 
of zero-measure circle-type MUPOs. 




FIG. 1: The drivebelt billiard is defined by the ratio of the small 
and big radii E (0, 1) and the arc of the larger of the two circu- 
lar arcs cf> S (tt/2, it) which in turns defines which MUPOs exist 
in the drivebelt's phase space. Left: Drivebelt with = 0.3 and 
<j> = 2.5. The period two diametrical MUPO (s, j) = (2. 1), the 
period three triangle MUPO (s,j) = (3,1) and the period four 
square MUPO (s,j) = (4,1) are shown. Right: These MUPOs 
can be oriented within a certain range depending on the value 
of (f>. A small hole is placed near the edge of the larger arc as to 
intersect the range of all 3 MUPOs present. 



Circular type MUPOs (defined later) in a sense gen- 
eralize bouncing ball orbits as they allow for peri- 
ods s > 2. They have also been observed in the 
mixed phase space annular |31j and mushroom bil- 
liards [32]. In fact, in both of these cases an infinite 
number of them was shown to exist [33] while care- 
ful investigations involving continued fractions and re- 
sults from number theory have revealed a MUPO-free 
set of mushrooms [32] which can also be extended to 
the annular billiard. In contrast, the drivebelt billiard 
has only a finite number of MUPOs controlled by the 
parameter (f>, hence motivating the term multiple in- 
termittency. That is, a finite number of regions of the 
billiard phase space individually display quasi-regular 
behavior for long periods of time. Interestingly and un- 
like bouncing ball orbits, circle type MUPOs may have 
large incidence angles with the boundary and therefore 
would be interesting to study in the context of micro- 
resonators [35) where total internal reflection will trap 
some of them according to Snell's law of refraction. Re- 



gardless of this, all aforementioned classical and quan- 
tum 'side-effects' due to bouncing ball orbits including 
nonstandard CLT of [T7] [T5] are expected to be inher- 
ited by the drivebelt. Moreover, an interesting duality 
exists between the stadium's bouncing ball orbits and 
the drivebelt's circular MUPOs as the roles of straight 
and curved boundary segments are reversed. Also, as 
we shall show, they exhibit very similar reflection laws 
at the end of each quasi-periodic intermittent sequence 
of reflections, a key property in our approach for cal- 
culating the asymptotic survival probability. 

In the next section we formally present our model by 
defining the relevant parameters and briefly discussing 
how MUPOs affect the open and closed billiard dy- 
namics. In section III, we numerically investigate the 
open drivebelt phase space and observe the multiple in- 
termittency manifested by the multiple sticky regions 
surrounding the system's MUPOs. In section IV we 
adapt and extend the approach of Ref [25] and obtain 
an exact to leading order in t analytic expression for 
the asymptotic survival probability function which we 
numerically verify and discuss different scaling limits. 
In section V we conclude with a short discussion. 



II. CIRCLE-TYPE MUPOS AND MULTIPLE 
INTERMITTENCY IN THE DRIVEBELT 
BILLIARD 

We define Birkhoff coordinates (z,sinf?) such that 
z E [0, \ dQ\) increases from zero anticlockwise from 
the rightmost point of the billiard as shown in Figure 
[T] and 9 e (— f , f ) is the angle of incidence. These 
coordinates are volume preserving under the billiard 
map and are also the ergodic equilibrium measure of 
the closed billiard and corresponding volume element 
dfi = (2\dQ\)^ 1 cos 6* dddz which we use as initial par- 
ticle ensemble in our numerical simulations since the 
hole will typically be of small but finite size. 

The drivebelt billiard is a non-uniformly hyperbolic 
system due to the existence of circle-type MUPOs in 
the larger of its two circular arcs. MUPOs in the drive- 
belt are a subset of periodic orbits of the circle billiard 
and hence have incidence angles given by 6 S j = f — 
where s > 2 and j > 1 are coprime integers describ- 
ing the period and rotation number respectively. The 
number of MUPOs in the drivebelt grows approxi- 
mately as ~ 3/(2(7r — 0) 2 ), hence we coin the term 
'multiple intermittency'. This is because a MUPO of 
period s will exist only if it satisfies — > 2(7T — <j>) and 
the number of integers coprime to s goes like ~ 6s /ir 2 . 
We therefore define the set containing all MUPOs 
for a given </> by 5 := {(2, 1), (s, j) 1 3 < s < 

1<J< L^J ,gcd(«,i) = l}. 

The mountain-plot in Figure [2j shows how MUPOs 
come into existence as cj> is increased from ^ to n from 
top to bottom. Also in Figure [2] and the right panel 
of Figure [l] we can see how the families of MUPOs 
define a range along which they can be oriented inside 
the larger of the two circular arcs. The example case 
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FIG. 2: As (j> S (7r/2,7r) is increased, more and more MUPOs 
come into existence. Only the first 7 MUPOs are shown here. 
The total range they cover also increases monotonically with 
4> and is shown on the horizontal axes. The range covered by 
MUPOs for the example case of </> = 2.5 is explicitly shown (also 
see right panel of Figure III. 



shown in the right panel of Figure [TJwith <f> = 2.5 is 
also highlighted in Figure [2] as a red horizontal line. 

The linear stability (monodromy matrix) eigenval- 
ues of MUPOs are both equal to one hence any small 
perturbation of a MFJPO will initially grow linearly 
(not exponentially) with time as the orbit precesses 
in the perturbed direction. Some examples are shown 
in Figure [l] MUPOs remain periodic after a small 
perturbation in z and hence form continuous families 
while perturbations in 9 will cause them to precess and 
exhibit long periods of quasi-regular behavior (approx- 
imately periodic), before leaving the vicinity of the 
periodic orbit and exploring the rest of the ergodic 
phase space uniformly and exponentially fast. This 
quasi-regular, non-dispersive behavior can cause the 
formation of scars in quantum billiards and is of par- 
ticular interest in optical microlasers where total inter- 
nal reflection can pick out some of the MUPOs while 
allowing others to cause directional emission |35j . 

Although MUPOs occupy zero volume in phase 
space and hence do not affect the overall ergodicity 
of the system, they govern long time statistical prop- 
erties of the closed system, such as the rate of mix- 
ing (the rate of the decay of correlations) C(t) ~ t^ 1 
[2"B] and the Poincare recurrence times distribution 
Q(t) ~ t~ 2 [36] which is also intimately related [37] 
with the long time survival probability P(t) ~ t~ x [35] 
of the open system. Furthermore, the exponents of 
these power-laws appear to be a universal fingerprint 
of non-uniform hyperbolicity and stickiness, at least 
for one and two dimensional area preserving maps with 
sharply divided phase space [3"5] . 

To this end, we note that the large number of consec- 
utive collisions that occur in a MUPO's close vicinity 
in phase space does not violate Kac's famous lemma 
[39j . That is, the mean recurrence time to a small box 
surrounding a MUPO equals one over the measure of 
the box. This is of course nothing more than an ar- 



tifact of ergodicity in that the long times spent by a 
trajectory stuck in the vicinity of the MUPOs are on 
average compensated by an apparent 'difficulty' in en- 
tering such sticky regions in the first place [2B]. Note 
however that due to the sticky dynamics, the second 
and higher moments of the recurrence times distribu- 
tion diverge [4*0] . 

The long periods of quasi-regular behavior and the 
ergodic compensation mechanism are particularly in- 
teresting to discuss in the multiple intermittency set- 
ting of the drivebelt billiard. This is because by defin- 
ing conditional return and transmission distributions, 
one may succeed in breaking the usual independence 
between events observed in chaotic systems. For ex- 
ample, consider small boxes A, B, etc. in phase space 
enclosing the MUPOs. Mixing of this billiard says that 
the probability of being in A at some time and B at 
a much later time approaches fx(A)fi(B). A priori, it 
does not say what happens at closely spaced times, 
and in particular whether the probability of reaching 
B rather than C immediately following A is propor- 
tional to their respective measures. If we require a 
given number n > 1 collisions within each region, the 
measures are reduced according to the collision rules 
where the MUPO reaches the straight segments, which 
leads to a dependence on the period of the orbit. Pre- 
liminary investigations (omitted) suggest that these 
conditional transmission probabilities are remarkably 
close to the default independence prediction, but de- 
viations must be evident at some level, in a way that 
reveals deeper structure of the dynamics of this, or 
indeed other, system with multiple intermittency We 
leave further investigation of this point to future work, 
but note that such a situation is particularly relevant 
to optical microlasers where MUPO resonances can in- 
teract via different mechanisms [551 FIl~rl4"5] and a min- 
imum number of collisions inside the cavity is needed 
to achieve high quality factors. In the next section we 
perform a numerical investigation of the open drive- 
belt phase space and identify the set of sticky orbits 
responsible for the power-laws discussed above. 



III. NUMERICAL INVESTIGATIONS OF 
DRIVEBELT PHASE SPACE 

We begin our investigations with some phase space 
images of a drivebelt billiard with R = 1, r = 0.3 and 
<j) = 2.5. Notice that = {(2, 1), (3, 1), (4, 1)}. Fig- 
ure [5] shows 10 5 randomly chosen initial conditions 
(ICs) (initially distributed on the billiard boundary 
according to the equilibrium density) which do not es- 
cape after time t = 2000 through a small hole of size 
e = 0.1 placed at z € (2.8, 2.9). Note that the hole is 
placed well inside the ergodic component of the phase 
space but is also away from any sticky regions. The 
hole can be clearly identified as the thin gray verti- 
cal strip on the right of the figure. We notice that 
long surviving initial conditions seem to populate the 
sticky regions surrounding the locations of the MU- 
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FIG. 3: Phase space plot of drivebelt billiard showing 10 5 
randomly chosen initial conditions which survive up to time 
t = 2000 with drivebelt parameters <j> = 2.5, R = 1, r = 0.3 
and a hole of size e = 0.1 positioned at z £ (2.8, 2.9). 



POs. We also notice that these ICs are supported on 
spike-like fractal-looking structures which correspond 
to the unstable manifold of the hole up to t = 2000. 
The coarse grained ICs surrounding these spikes cor- 
respond to trajectories that move from one MUPO to 
another but do not intersect the hole up to t = 2000. 
Moreover, these long surviving heteroclinic orbits con- 
tribute to the asymptotic power-law decay even though 
they exhibit stable and unstable directions in phase 
space. 



in and take a closer look (see Figure pj. The magni- 
fied phase space corresponds to z € (2^45, 2.5) which 
is right where the larger circular arc meets with the 
straight. The spike-like fractal- looking structure ob- 
served previously is retained for negative 9 but is lin- 
ear for positive 9. This pattern is preserved in all 
MUPO vicinities (not shown) suggesting that for long 
but fixed survival times t, small and positive 9 pertur- 
bations of MUPOs have a linear dependence with dis- 
tance to the hole while negative 9 perturbations have 
a more interesting dependence and can also sustain 
larger perturbation magnitudes. In the next section 
we clarify these two types of behaviors at long times 
by classifying them as ICs precessing towards/away 
from the hole. 
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FIG. 5: Magnification of phase space plot in Figure [4] at 2 £ 
(2.45,2.5) and sin0 6 (-0.0004,0.0004) showing how the long 
surviving t > 2000 initial conditions populate the area near the 
period two diametrical MUPO. 
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FIG. 4: Phase space plot of drivebelt billiard showing 10 5 
randomly chosen initial conditions which survive up to time 
t = 2000 with drivebelt parameters <j> = 2.5, R = 1, r = 0.3 
and a hole of size e = 0.1 positioned at 2 6 (2.3, 2.4). 



Figure [4] shows 10 5 randomly chosen initial condi- 
tions which do not escape after time t = 2000 through 
a small hole of size e = 0.1 intentionally placed at 
z e (2.3, 2.4) as to overlap part of all three MUPO 
ranges. Note that the number of initial conditions 
plotted in Figures [3] and [4] are the same while now 
the hole clearly repeats itself in all sticky regions. The 
density of long surviving orbits in regions away from 
MUPOs appears to have decreased while the spike like 
structure has shrank dramatically. We therefore zoom 



IV. ASYMPTOTIC SURVIVAL 
PROBABILITY 

Assuming that for large enough times the long sur- 
viving initial conditions are dense only in the close 
vicinity of MUPOs, then a hole which overlaps part of 
all MUPO ranges will separate the long surviving ICs 
into two simple families: 1) orbits initially precess- 
ing towards the hole, and 2) orbits initially precessing 
away from the hole. We clarify this classification as 
follows. A small perturbation rji <C 1 in the angle of 
incidence 9 s j of a MUPO will cause the orbit to pre- 
cess in a quasi-periodic fashion. If the orbit is precess- 
ing towards the hole and rji is small enough then the 
orbit may survive for an unbounded amount of time 
and hence will contribute to the asymptotic power-law 
~ t^ 1 tail of P(t). In fact such orbits correspond to the 
ICs with 9 > shown in Figure [5] and are bounded lin- 
early from above as the precessing velocity is propor- 
tional to the perturbation strength. If the orbit is pre- 
cessing away from the hole and the % is small enough, 
the orbit will experience a nonlinear collision processes 
which reverses the orbit's precession direction and is 
of the type C ... CSC S - X (S)C ...C until escape. We 
have used here a symbolic dynamics where the sym- 
bols S and C correspond to collisions on straight and 
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curved boundary segments respectively. The s in the 
superscript corresponds to the period of the nearest 
MUPO in phase space. Also the S in brackets, as we 
shall see later only occurs when the first straight seg- 
ment collision does not reverse the precessing orbit's 
direction. After the precessing direction is reversed, 
the orbit will approach the hole with a final angle of 
incidence which is not far from 9 s j and may again 
have unbounded survival time. Such ICs as we shall 
soon see constitute the spike-like structure observed 
for 9 < in Figure [5j With this classification, we 
expect that 



lim tP(t) = lim i(Pi(i) 

t-^-oo t—>oo 



(2) 



where Pi (t) and P2 (t) correspond to the leading order 
contributions of the two types of ICs. 



A. Precessing towards the hole 

Consider a small hole of size e situated at z £ 
(h~,h + ) C (0, </>), such that it overlaps all MUPO 
ranges and h + = h~ + e. An example is shown in 
the right panel of Figure [T] We perform a change 
of variables and consider ICs (4>i,9i) such that <fo £ 
(0, (j) — h + ) is the angular distance of the IC from the 
nearest edge of the hole in the direction of precession 
and 9i — 9 s j + r\i is the angle of incidence such that 
the perturbation rji <C 1 causes the orbit to precess 
towards the hole. Since the precessing angular veloc- 
ity is proportional to the perturbation strength rji the 
condition 



\m\ < 

2s 



(3) 



guarantees that the IC does not 'jump over' the hole. 
Hence such an orbit will escape in time t given by 



t( 



i) 



2^s 



R4>i cos 9 s j 



Vi 



s2R cos 9i 

0(1). 



(4) 



Notice that we have obtained the escape time in a con- 
tinuous time frame from discrete iterations of the bil- 
liard map. This is allowed here since the time between 
collisions is fixed. From the above information we can 
obtain a lower bound on the relevant time scales for 
which our results are valid 



\t]i\ < mm 



e R4>i cos 9 s j 



2s' 



and hence when 



t > 



4vri? 



(5) 



(G) 



since sup(0i) = 2ir/s. 

We can now integrate the area described by the 
above inequalities with respect to the equilibrium 



measure (2\dQ\)- 1 dzdsm6 = (2\dQ\)~ 1 R(cos9 sJ 
rji sin 9 s j)d(j)idr]i as follows: 



i?(cos s j — r/i sin 9 S 



2\dQ\ 



i? 2 ^(3t - Rfa sind sJ ) cos 2 0. 
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(7) 



where c/H(s) is the angular distance available within 
the range of the MUPO on the corresponding side of 
the hole. Similarly we obtain the asymptotic contri- 
bution from ICs precessing towards the hole but from 
the other side of the hole 



_ i^cos 2 ^ (1 



(8) 



and the total contribution to the asymptotic survival 
probability is 



(9) 



where the p — 1 if s = 2 and p — 2 otherwise as to ac- 
count for the vertical symmetry of the phase space and 
s accounts for the multiplicity of the MUPOs as seen 
in Figures [3] and [4j We next turn to the contribution 
P2W due to ICs precessing away from the hole. 



B. Precessing away from the hole 

We now consider ICs (4>i,9i) such that the pertur- 
bation r/i <^ I causes the orbit to initially precess away 
from the hole. Here, we define & G (0, cj)—h + ) to be the 
angular distance of the IC from the edge of the circu- 
lar arc at z = (j). The precessing orbit will first collide 
with the straight segment after sn £ N + consecutive 
C-type collisions with the boundary during which time 
9i is constant, n here is the number of MUPO periods 
the IC will complete just after traversing the available 
curved distance to z = <j). The first «S-type collision will 
occur at z — 4>+di where di = 2\rji\sn— 4>i~\-0(rjf) > 0. 
Since the collision occurs on a straight segment of the 
billiard's boundary, we may 'unfold' the billiard using 
the image reconstruction trick [32] (see Figure [6]). The 
resulting angle of incidence with the curved segment is 
to leading order 9f = 9 t — 2d\ = 9 s j + r/f. Hence, we 
conclude that the nonlinear collision process described 
above causes an orbit initially perturbed by rji to have 
a final incidence angle 9f which is r/f away from the 
original MUPO angle Sij such that 2i = 1 - Ml., if 
2 1 rjf I s > g?i then the orbit will start to precess towards 
the hole and will eventually escape without jumping 
over the hole if 2|?y/|s < e. The later is guaranteed for 
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FIG. 6: A slightly perturbed in angle MUPO with initial angle 
of incidence 8i = 8 a j + r\i is 'unfolded' when it collides on a 
straight boundary segment giving a new angle of 8 f = 8i — 2di = 
s ,j +Vf- 



t large enough. The symbolic dynamics characteriz- 
ing the trajectory of such an IC is therefore given by 
C . . . CSC . . . C. If however 2|r//|s < d\, then a second 
5-type collision will occur such that the final angle is 
now Of = 9 s j +rf f = 9 sJ + 4di(l - 2s) + rji(As - 1). 
The symbolic dynamics sequence of such a collision 
process is hence given by C . . . CSC S ~ 1 SC . . . C . The 
appropriate scaling which linearizes \rjf/rji\ is given by 
A = • Using this we arrive at a small perturbation 



collision rule given by 




FIG. 7: The dependence of | on A for the period two diamet- 
rical MUPO (s = 2), where ?7i <gC 1 and A = as described 



in equation (| 1 Op 



2>7 4 s 



We now formulate the time to escape function: 



Rcf>i cos @s,j 



s2i?cos 6 



2\nj\ 8 
Rfa cos 9 s j 



s2Rcos6f + S 



(11) 



to leading order where 5 = s2R cos 9 s j +0(r]i) and we 
have assumed that both rji and r/f are small. Notice 
that again we have obtained the escape time in a con- 
tinuous time frame from discrete iterations of the bil- 
liard map. This is allowed here since the time between 
collisions is fixed and deviates only slightly during the 
nonlinear collision process. We drop the modulus signs 
from here on, neglect higher order terms in rji and r/f, 
and rearrange to obtain: 



= Rr/ftfii cos I 



Rr]i(f>i cos 9 s j - tr)iT)f, (12) 



which is a conic section since it is quadratic in both fa 
and rji due to the definition of d\. Substituting rjf into 
(12 1 and writing everything in terms of fa and 7jj, we 



obtain two hyperbolas for each value of n. Each pair of 
hyperbolas encloses an area which stretches and tilts 
in a non-overlapping fashion as n is increased. Figure 
[8] shows how this approximation reproduces the area 
occupied by the long surviving orbits in the vicinity of 
the diametrical period two MUPO. 



-16s 2 A + s(4 + 8A) - 1, 
4As - 1, 



0<A<4^T 
< A < 1 



4s-l 



/ 
/ 

(10) 



= 2 
= 1. 



where / = 1, 2 is the corresponding number of curved 
collisions during the nonlinear collision process. This 
is plotted in Figure [7] for the period two diametrical 
orbit. Notice that when \rjf/rji\ > 1, the return pro- 
cessing speed of the IC is increased and occurs with 
higher probability than a decrease. Also, a double 
(/ = 2) collision scenario is much less probable than a 
single one (/ = 1). 




FIG. 8: The first 100 spikes produced by the hyperbolas of 
equation i ] 1 2[ ) for the same parameter values as in Figure [5] 
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Noticing that the hyperbolas defined by (12 1 look 
like straight lines in Figure |SJ we use the intersec- 
tion coordinates at the top and bottom of each spike 
to form lines which maintain the general structure of 
the enclosed surviving set. Note that the error in this 
approximation is only of order 0(t~ 2 ). This allows for 
easy integration and summation over the area of each 
spike in order to obtain Pz(t) to leading order. How- 
ever, care must be taken at large values of n since the 
spikes will tilt into the hole and hence modifications 
are necessary. For a step- by-step explanation see [25] . 

Finally, adding all MUPO contributions and ex- 
panding for large t we obtain 



p 2 (t)= J2 s p 

(4a - 1) 



R 2 {<t>\ + cos 2 sJ 
4|0Q|i 



(13) 



x fl+ ln(4s - 1) 

L 2a(2a -1) K . 



where p was defined in the previous subsection and 
<j)\ and 4>2 are the allowed ranges on either side of the 
hole and hence implicitly depend on s. Note that this 



equation is valid for times t > t = 



8(s0-7r) 



which is 



larger than ^ and can be derived in a similar way. 



C. Asymptotic P(t) and numerical simulations 




FIG. 9: Log-linear plot of the Survival probability P(t) as a 
function of time t for the open drivebelt with <j> = 2.5, R = 1, 
r = 0.3 and a hole of size e = 0.1 positioned at z £ (2.3, 2.4) as 
can be seen in the inset. The blue curve is given by the analytic 
prediction of equation |15| while the empty black circles are 
from a numerical simulation consisting of 10 s initial conditions. 
The inset is a plot of tP(t) showing the agreement with the 
analytic expressions for the constant C. 



In the previous subsections we first classified long 
surviving ICs into two simple families: 1) orbits ini- 
tially precessing towards the hole, and 2) orbits ini- 
tially precessing away from the hole. This simplifica- 
tion was achieved by placing the hole as to overlap all 
MUPO ranges. We then derived conditions defining 
the sets of the ICs surviving up to time t for each clas- 
sification and integrated their volume in phase space. 



Adding the two contributions and taking the limit we 
obtain for the survival probability of the open drivebelt 
with initial conditions chosen according to the invari- 
ant measure on the boundary 



C = lim tP(t) = ^2 sp 



R 2 {(t)\ + (/>!) cos 2 e Stj 



x 2 + 



(4* - 1) 
2s(2s - 1 



4|5Q| 
ln(4s - 1) 



(14) 



We now confirm our theoretical predictions through 
numerical simulations. These are shown in Figure 
[9] indicating an excellent agreement with the theory. 
Other simulations with different <f> and e have also been 
performed but are not included here. Moreover, the 
numerical simulations suggest that the survival prob- 
ability function exhibits a coexistence of exponential 
and power law decay as in typical weakly chaotic sys- 
tems [44] given by: 



P(t) = 



irregular, 



for t <t 
for t > i , 



(15) 



where we have neglected terms of order ~ t~ 2 . The 'ir- 
regular' short-time behavior shows small fluctuations 
which are a result of geometry dependent short orbits. 
These become less important if the hole is small. Also, 
the coefficient of the exponential term is approximately 
1 since for small holes, mixing causes the system to 
forget its initial state and therefore the probability de- 
cays as a Poisson process. Finally, equation (15) sup- 



ports the arguments of [12 in that the chaotic saddle 
splits into a hyperbolic (exponential) and nonhyper- 
bolic (power-law) part which give separate contribu- 
tions to P(t). 



D. Scaling 

The natural scaling limit for open billiards (or in- 
deed other dynamical systems) is when r = et is kept 
constant while e — > and t — > oo. Substituting this 
into ([!]) and (151 above we get that 



lim P (We) 

e->0 



(16) 



since C/t = eC/r with C constant vanishes in this 
limit. Hence, the open drivebelt billiard (and also 
the straight stadium which has a similar form of the 
survival probability) satisfies an exponential scaling 
limit in common with more chaotic systems of ([!]). In 
other words such scaling approaches "miss" the asymp- 
totic tails and cannot discriminate between strongly 
chaotic and weakly chaotic systems. It is worth men- 
tioning that this exponential scaling limit is also in- 
timately connected with limiting distributions of re- 
currence times [45j [46] . An exception to exponential 
behaviour is the completely regular circular billiard 
which shows (numerically) a nontrivial scaling law [3] . 
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V. DISCUSSION AND CONCLUSIONS 

We have considered the open drivebelt billiard and 
generalized the approach of Ref.[25] in order to ob- 
tain exact to leading order expressions for the asymp- 
totic survival probability P(t) ~ C/t. We have con- 
firmed our predictions through numerical simulations 
and have obtained a full description of P(t) to lead- 
ing order. The generalization of our previous results is 
based on a small perturbation collision rule derived 
at the vicinity of smoothly connected circular and 
straight boundary segments. Therefore, our approach 
is now complete and can further be applied to a variety 
of smooth billiards constructed of conic components 
and straight lines (for example elliptical stadia). 



The quantitative result offered by equation ( 14 1 and 



its derivation goes beyond previous investigations of 
open billiards which were only of qualitative nature in 
that it provides answers to question like where to place 
the hole achieve maximum / minimum escape |47j . This 
can be done by minimizing or maximizing C as a func- 
tion of hole position. For the drivebelt, we find that 
the maximum C is attained when the hole is placed 
right at the edge of the arc i.e. at z € ((/> — e, <j>). The 
minimum is a more complicated function of <f>. 

Our approach for obtaining the constant C is how- 
ever restricted to the case of a hole positioned as to in- 
tersect all MUPO ranges and thus splitting the chaotic 
saddle into non-hyperbolic and hyperbolic parts such 
that the former occupies a simple region in phase space 
which can be isolated and integrated over. The case 
of a hole not intersecting all MUPO ranges remains 
an open problem although the power-law asymptotic 
decay ~ t^ 1 is expected to persist. Interestingly, the 
multiple intermittency of the drivebelt will cause mul- 
tiple escape paths to exist. For instance, a sticky 



orbit may exit and then enter into a different sticky 
mode with some transfer probability before escape oc- 
curs. Such scenarios motivate investigations of escape 
through multiple holes where different mechanisms can 
be utilized to produce for example asymmetric trans- 
port as observed in the stadium billiard [48j . Similarly, 
one may consider how an increased number of holes in- 
tersecting different MUPOs can affect escape through 
shadowing effects 03] ■ Such considerations are very 
relevant in the field of controlling chaos since the es- 
cape route through one of the holes may be consider- 
ably reduced by other holes. 

The results on the asymptotic survival probability 
in the open drivebelt and the improved understanding 
of stickiness due to MUPOs presented in this paper are 
hoped to shed new light on the otherwise not so well 
understood quantum analogue of stickiness (scarring 
and anti-scarring). Similarly, we hope to extend our re- 
sults in a semiclassical regime thus offering important 
ramifications in the context of quantum chaos with 
experimental applications such as electronic transport 
through open ballistic nano-structures (quantum dots) 
[50] . Finally, by defining conditional recurrence time 
distributions we expect to be able to manipulate and 
adapt the multiple intermittency in the drivebelt as 
to investigate physically interesting scenarios such as 
achieving high quality factors and directional emission 
from dielectric cavities @]- 
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